setwd("C:/Users/cliente/Dropbox/UFMT/Disciplinas/Veterinaria/experimental/dados")

setwd("D:/experimental/dados")
require(xlsx)
dados=read.xlsx("insulina.xls",1)
attach(dados)

####Gerar a análise de variância
model=aov(VR~trat)
anova(model)


##########################################################################
#####################Comparações multiplas################################
##########################################################################
require(multcomp)
##Teste t
contrast<-rbind(c(1,-0.5,-0.5),c(1,0,-1))
contrast
summary(glht(model,mcp(trat=contrast)))


#Teste de Tukey
require(agricolae)
HSD.test(model,"trat", group=TRUE)

#Teste SNK
gl=df.residual(model)
QME=deviance(model)/gl
SNK.test(VR,trat,gl,QME, group=TRUE)

#Teste de Duncan
duncan.test(VR,trat,gl,QME, group=TRUE)


#Teste de Scheffé
zz=anova(model)
Fc=zz$F[1]
scheffe.test(VR, trat, gl,QME, Fc, alpha = 0.05, group=TRUE)

#Teste de Dunnett
require(multcomp)
grupos=tapply(VR,trat,length)
grupos
contrast=contrMat(grupos, base=2,type ="Dunnett") ##base é igual coluna do controle
constrat
summary(glht(model,mcp(trat=contrast)))

require(ScottKnott)
teste=SK(trat, VR, model,which='trat', sig.level=0.05)
summary(teste)